Monte Carlo simulation of single spin asymmetries in pion-proton collisions 
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We present Monte Carlo simulations of both the Sivers and the Boer-Mulders effects in the po- 
larized Drell-Yan ir — > ^ + f/T X process at the center-of-mass energy y^i ~ 14 GeV reachable 
at COMPASS with pion beams of energy 100 GeV. For the Sivers effect, we adopt two different 
| parametrizations for the Sivers function to explore the statistical accuracy required to extract un- 

ambiguous information on this parton density. In particular, we verify the possibility of checking 
its predicted sign change between Semi-Inclusive Deep-Inelastic Scattering (SIDIS) and Drell-Yan 
processes, a crucial test of nonperturbative QCD. For the Boer-Mulders effect, because of the lack 
, of parametrizations we can make only guesses. The goal is to explore the possibility of extracting 

• information on the transversity distribution, the missing piece necessary to complete the knowledge 

Pjr^ ' of the nucleon spin structure at leading twist, and the Boer-Mulders function, which is related to 

the long-standing problem of the violation of the Lam- Tung sum rule in the unpolarized Drell-Yan 
cross section. 

PACS numbers: 13.75.Gx,13.88+e,13.85.Qk 
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^ ! I. INTRODUCTION 

O - 

\Q . The recent measurement of Single-Spin Asymmetries (SSA) in semi-inclusive lp' — > I'irX Deep-Inelastic Scattering 
(SIDIS) on transversely polarized hadronic targets 0- d []| ■ has renewed the interest in the problem of describing the 
spin structure of hadrons within Quantum Chromotypes (QCD) ||, and has stimulated since then a large pro- 
O ' duction of phenomenological and theoretical papers. Experimental evidence of large SSA in hadron-hadron collisions 
was well known since many years 0, , but it has never been consistently explained in the context of perturbative 
QCD in the collinear massless approximation The idea of going beyond the collinear approximation opened new 
perspectives about the possibility of explaining these SSA in terms of intrinsic transverse motion of partons inside 
hadrons, and of correlations between such intrinsic transverse momenta and transverse spin degrees of freedom. The 
most popular examples are the Sivers [fj and the Collins [Tol ] effects. In the former case, an asymmetric azimuthal 
distribution of detected hadrons (with respect to the normal to the production plane) is obtained from the nonper- 
turbative correlation p T x P ■ S T , where p T is the intrinsic transverse momentum of an unpolarized parton inside 
a target hadron with momentum P and transverse polarization S T . In the latter case, the asymmetry is obtained 
from the correlation k x P} lT ■ s T , where a parton with momentum k and transverse polarization s T fragments into 
an unpolarized hadron with transverse momentum P} lT . In both cases, the sizes of the effects are represented by new 
Transverse-Momentum Dependent (TMD) partonic functions, the socalled Sivers and Collins functions, respectively. 
However, SSA data in hadronic collisions have been collected so far typically for semi-inclusive 

pp m h^X 

processes, where the factorization proof is complicated by higher-twist correlators and the power-suppressed 
asymmetry can be produced by several (overlapping) mechanisms. On the contrary, the theoretical situation of the 
SIDIS measurements is more transparent. On the basis of a suitable factorization theorem [T3. fT3| . the cross section 
at leading twist contains convolutions involving separately the Sivers and Collins functions with different azimuthal 
dependences, sin(0 — <j>s) and sin(^ + 4>s)i respectively, where <fi, <ps, are the azimuthal angles of the produced hadron 
and of the target polarization with respect to the axis defined by the virtual photon ^4[. According to the extracted 
azimuthal dependence, the measured SSA can then be clearly related to one effect or the other 0, Q . 

Similarly, in the Drell-Yan process U\Ii\ — > l + l~ X the cross section displays at leading twist two terms weighted by 
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sin(0— 4>s) and sin(0+0s), where now 0, 0s, are the azimuthal orientations of the final lepton plane and of the hadron 
polarization with respect to the reaction plane |l5j . Adopting the notations recommended in Ref . [T^ | , the first one 
involves the convolution of the Sivers function with the standard unpolarized parton distribution fi. The second 
one involves the transversity distribution hi and the Boer-Mulders function h^, a TMD distribution which is most 
likely responsible for the violation of the Lam- Tung sum rule in the corresponding anomalous cos 2<p asymmetry of 
the unpolarized Drell-Yan cross section |Tfij . Hence, a simultaneous measurement of unpolarized and single-polarized 
Drell-Yan cross sections would allow to extract all the unknowns from data [TtI fl8| . Both hi and hj- describe the 
distribution of transversely polarized partons; but the former applies to transversely polarized parent hadrons, while 
the latter to unpolarized ones. On an equal footing, and fi describe distributions of unpolarized partons. The 
correlation between p T and S T inside f^ T is possible only for a nonvanishing orbital angular momentum of partons. 
Then, extraction of Sivers function from SIDIS and Drell-Yan data would allow to study the orbital motion and the 
spatial distribution of hidden confined partons [l9| , as well as to test its peculiar universality property |20j| . 

In a series of previous papers, we performed numerical simulations of single-polarized Drell-Yan SSA for the pp^ — > 
[i + H~X [2l| and pp^ — > ^ [>T X 0] processes. With proton beams, we considered collisions at ^/s — 200 GeV in 
the kinematic conditions for the foreseen upgrade of RHIC (RHIC II) . Even if in pp collisions the nonvalence partonic 
contribution to the elementary annihilation is unavoidable (leading, in principle, to lower counting rates), still the 
kinematics selects a portion of phase space that emphasizes this contribution. The net result is that with a reasonable 
sample of Drell-Yan events the statistical accuracy allows to unambiguously extract the Sivers function from the 
corresponding sin(0 — cf>s) asymmetry, as well as to clearly test its predicted sign change with respect to the SIDIS 
asymmetry [21|. In pp collisions, the cross section is dominated by the valence contribution to the annihilation of 
a parton (from p) and an antiparton (from p); hence, in general it is not suppressed as in the previous case (for a 
quantitative check in our Monte Carlo, see Sec. IVB of Ref. jlj). In Ref. we selected antiproton beams of 15 
GeV, as they could be produced at the High Energy Storage Ring (HESR) at GSI [22ll2^ |. and we simulated collisions 
at y/s ~ 14 GeV in the socalled asymmetric collider mode. The goal was to explore the minimal conditions required 
for an unambiguous extraction of hi and h^ from a combined analysis of the sin(</> + (j>$) and cos 20 asymmetries in 
the full (unpolarized + polarized) cross section. 

Here, we will reconsider the same scenarios but for the 7r ± pT _> ^+ jU - x process at the same i/s ~ 14 GeV that can 
be reached at COMPASS with pion beams of energy 100 GeV. As for p beams, the elementary mechanism is dominated 
by the annihilation between valence partons (from p) and valence antipartons (from it) . Indeed, a large Sivers effect 
was predicted in this context by usig the same Sivers function fitted to the measured sm{4> — <ps) asymmetry in 
SIDIS 24] . Taking advantage on the high statistics reachable with pions, in our Monte Carlo we simulate both 
sin(0 ± (f>s) SSA in the Drell-Yan cross section. For the Sivers effect we use two parametrizations of fx T : the one of 
Ref. HH, which was deduced by fitting the recent HERMES data for the sin(0-0 s ) SSA 0; the one of Ref. [2l|, which 
is constrained by the recent RHIC data for the pp^ — > irX process at higher energy [2£j. For the Boer-Mulders effect, 
since there is no such abundance of data and fits, we follow Ref. [15J to constrain h± by the azimuthal asymmetry of 
the corresponding unpolarized Drell-Yan cross section (see also Ref. [27I for a similar analysis). Then, we insert, as 
we did in Ref. |17| . very different input test functions for hi in order to explore the sensitivity of the simulated SSA 
within the statistical accuracy. 

In Sec. [n] we review the formalism and the details of the numerical simulation. In Sec. IIIII we present and discuss 
our results. Finally, in Sec. II VI we draw some conclusions. 




FIG. 1: The Collins-Soper frame. 



II. GENERAL FRAMEWORK FOR THE NUMERICAL SIMULATION 



In a Drell-Yan process, an antilepton-lepton pair with individual momenta fci and ki is produced from the collision 
of two hadrons with momentum Pi, mass Mi, and spin Si, with i = 1,2. The center-of-mass (cm.) square energy 
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available is s = (Pi + P2) 2 and the invariant mass of the final lepton pair is given by the time-like momentum transfer 
q 2 = M 2 = (k\ + k 2 ) 2 - If M 2 ,s — > 00, while keeping the ratio < t = M 2 / s < 1 limited, a factorization theorem 
can be proven [28j ensuring that the elementary mechanism proceeds through the annihilation of a parton and an 
antiparton with momenta pi and P2, respectively, into a virtual photon with time-like momentum q 2 . If P± and 
P 2 are the dominant light-cone components of hadron momenta in this regime, then the partons are approximately 
collinear with the parent hadrons and carry the light-cone momentum fractions < x\ = Pi jP^ , X2 = P2 /P2 ^ 1; 



with q 



+ _ 



Pi , q 



p 2 by momentum conservation |15l |. The transverse components pi T of pi with respect to the 



direction defined by Pi(i = 1,2), are constrained again by the momentum conservation q T = p± T + P2t, where q T 
is the transverse momentum of the final lepton pair. If q T 7^ the annihilation direction is not known. Hence, it is 
convenient to select the socalled Collins-Soper frame [2!j described in Fig. ^ The final lepton pair is detected in the 
solid angle (9, 4>), where, in particular, (f> (and all other azimuthal angles) is measured in a plane perpendicular to the 
indicated lepton plane but containing h = Qt/|Qt|- 

The full expression of the leading- twist differential cross section for the ffiifj ~~ * l + l~X process can be written 
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where a is the fine structure constant, d£l = sin 9d9d(f>, e/ is the charge of the parton with flavor /, 
angle of the transverse spin of hadron i, and 



(1) 



!>Si is the azimuthal 
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The TMD functions // (H), h x (H), describe the distributions of unpolarized and transversely polarized partons in 
an unpolarized hadron H, respectively, while f^/(H^) and the pair h{(H^),hfJ(H^), have a similar interpretation 
but for transversely polarized hadrons W. The convolutions are defined as 



DF((Hi 



<df*(hP) 



dpi T dp 2T Hpit + P2t - q T ) DFi(xi,pi T ;f/Hi 



lDF2(x2,P2Tj/H^) + (f 



D 



(3) 

In previous papers, we made numerical simulations of the SSA generated in Eq. Q by the azimuthal dependences 
cos 2(j) and sin(0+0s 2 ) for antiproton beams H 1 = p by the sin(0— 4*s 2 ) dependence for proton beams H 1 = p [2ll |. 
as well as for double-polarized Drell-Yan processes with Pl\ = H\ = p^ (3(| . A combined measurement of these SSA 
allows to completely determine the intertwined unknown transversity hi and Boer-Mulders function h±, and the 
Sivers function f^. The Monte Carlo simulation was performed for high-energy proton beams (y/s = 200 GeV) in the 
conditions of the foreseen upgrade of RHIC (RHIC II), and for antiproton beams of 15 GeV as they could be produced 
at HESR-GSI. In the latter case, several scenarios were explored for 5 < < 14 GeV and 1.5 < M < 2.5, 4 < M < 9 
GeV, in order to avoid overlaps with the strange, charm, and bottom quarkonia [where the elementary annihilation 
does not necessarily proceed through a simple intermediate virtual photon, as it is assumed in Eq. Q]. Here, we 
reconsider the sin(</> — 4>s 2 ) an d sin(0 + (f>s 2 ) asymmetries by using pion beams of 100 GeV as they can be produced 
at COMPASS, in the fixed target mode such as to reach the same maximum cm. energy considered at HESR-GSI, 
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namely y/s ~ 14 GeV. Most of the technical details of the simulation are mutuated from our previous works; hence, 
we will heavily refer to Refs. 0, 0] in the following. 

The Monte Carlo events have been generated by the following cross section fl7| : 



da 



d n dxidX2dqT =K- s \r( qT , X u X2 ,M)f ± vfa^WbK) . ( 4 ) 

where the event distribution is driven by the elementary unpolarized annihilation, whose transition amplitude T has 
been highlighted. In Eq. I|T[1. we assume a factorized transverse-momentum dependence in each parton distribution 
such as to break the convolution J 7 , leading to 

\T\ 2 ^A(q T ,x 1 ,x 2 ,M)F{x 1 ,x 2 ), (5) 

where q T = \q T \- The function A is parametrized and normalized as in Ref. |3l| . where high-energy Drell-Yan n — p 
collisions were considered. The average transverse momentum turns out to be (q T ) > 1 GeV/c (see also the more recent 
Ref. [22J), which effectively reproduces the influence of sizable QCD corrections beyond the parton model picture of 
Eq. (Q. It is well known [33| that such corrections induce also large K factors and an M scale dependence in parton 
distributions, determining their evolution. As in our previous works p"?ll2ll,l30l |. we conventionally assume in Eq. Q 
that K = 2.5, but we stress that in an azimuthal asymmetry the corrections to the cross sections in the numerator 
and in the denominator should compensate each other, as it turns out to actually happen at RHIC cm. square 
energies [34j . Since the range of M values here explored is close to the one of Ref. [3l| , where the parametrization 
of A, F, and c,; in Eq. (0J was deduced assuming Af-independent parton distributions, we keep our same previous 
approach [l3)|2llll3 an d use 

2 

F{xu X2) = I2V E 4 /i>i;/M) /tea; f/H 2 ) + (f <- /) , (6) 

where the unpolarized distribution f( (x) for various flavors / = u,d, s, is taken again from Ref. |3l| . 

The whole solid angle (8, <fr) of the final lepton pair in the Collins-Soper frame is randomly distributed in each 
variable. The explicit form for sorting it in the Monte-Carlo is 0> 

Y / c l {q T ,x 1 ,x 2 )S l (9,cj ) ,<j ) s 2 ) = l + cos 2 fl+ S in 2 g cos 20 

1=1 

+ \S 2 T\c4(q T ,xi, X2 ) S4{9, 0, cf>s 2 ) ■ (7) 

If quarks were massless, the virtual photon would be only transversely polarized and the angular dependence would 
be described by the functions c% = Si = 1 and c 2 — 1, S 2 = cos 2 6. Violations of such azimuthal symmetry induced by 
the function 03 = ^ are due to the longitudinal polarization of the virtual photon and to the fact that quarks have an 
intrinsic transverse momentum distribution, leading to the explicit violation of the socalled Lam- Tung sum rule |3l| . 
QCD corrections influence v 7 which in principle depends also on M 2 [3ll |. Azimuthal cos 20 asymmetries induced 
by v were simulated in Ref. using the simple parametrization of Ref. and testing it against the previous 
measurement of Ref. [3l] . 

If we consider the Sivers effect in Eq. Q , the last term in Eq. becomes 

S 4 (0,0,0&) = (l + cos 2 #) sin(0-0 S2 ) (8) 

and the corresponding coefficient C4 reads 

f((x!,pi T ) fi T f (x 2 ,p 2T ) 



c±(q T ,xi,x 2 ) 



hp 2l 



M 2 



Ysf^fF //Oi,Pit) f{{x 2 ,p 2T ) 



(9) 



where the complete dependence of the involved TMD parton distributions has been made explicit. 
Viceversa, if we consider the Boer-Mulders effect in Eq. Q the last term in Eq. J7J) becomes 

5 4 (0, 0, 0s 2 ) - sin 2 9 sin(0 + Sa ) (10) 
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and the corresponding coefficient C4 reads 



04(^,0;!, x 2 ) 
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In the following, we will discuss different inputs for the x and p T dependence of these distributions which allow to 
calculate the convolutions and determine C4. In any case, following Refs. 0,0, [HI, the general strategy is to divide 
the event sample in two groups, one for positive values "U" of S4 in Eq. (JSJ) or i|l(J|) . and another one for negative 
values " D" , then taking the ratio (U — D)/(U + D). Data are accumulated only in the x 2 bin, i.e. they are summed 
upon xx,0, and q T . Statistical errors for the spin asymmetry (U — D)/(U + D) are obtained by making 10 independent 
repetitions of the simulation for each individual case, and then calculating for each x 2 bin the average asymmetry 
value and the variance. We checked that 10 repetitions are a reasonable threshold to have stable numbers, since the 
results do not change significantly when increasing the number of repetitions beyond 6. 



A. The Sivers effect 



Recently, the HERMES collaboration released new SSA data for the SIDIS process on transversely polarized pro- 
tons 0, which substantially increase the precision of the previous data set 0. As a consequence, different parametriza- 
tions of the Sivers function f^, have been extracted from this data set and found compatible also with the recent 
COMPASS data 3 (for a useful comparison among the various approaches see Ref. [35J). Following Ref. |2l|, we first 
simulate the Sivers effect using the parametrization of Ref. [25| , 



Sxt{x,Pt) = 
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(12) 



where M 2 is the mass of the polarized proton, p T = \p T \, and (pj) = 0.25 (GeV/c) 2 is deduced by assuming a Gaussian 
ansatz for the p T dependence of /1 in order to reproduce the azimuthal angular dependence of the SIDIS unpolarized 
cross section (Cahn effect). Flavor-dependent normalization and parameters in the x dependence are fitted to SIDIS 
SSA data neglecting the (small) contribution of antiquarks. The resulting parameters Mq and Nf, a,f,bf, with / = u,d, 
are listed in Tab. [I] The sometimes poor resolution of the fit forced us to select only the central values in order to 
produce meaningful numerical simulations. 

Following the steps described in Sec. III-l of Ref. j^, in particular the predicted sign change of f± T when going 
from SIDIS to Drell-Yan, we insert the opposite of Eq. (|12f> into Eq. @ and simplify it down to 



c 4 



4Mp q T 
& + 4M a 



(a u + b u ) 



a„+b„ 



a u ib u 



■x? {l-x 2 ) K + N, 
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(13) 



TABLE I: Parameters for the Sivers distribution from Ref. 25] 



quark up 




quark down 




N u 


0.32 ±0.11 


N d 


-1.0 ±0.12 


a v 


0.29 ±0.35 


da 


1.16 ±0.47 


b u 


0.53 ±3.58 


b d 


3.77 ±2.59 


Ml 


0.32 ± 0.25 (GeV/c) 2 







As an alternative choice, we adopt the new parametrization described in Ref. [2l|. It is inspired to the one of 
Ref. [13, where the transverse momentum of the detected pion in the SIDIS process was assumed to come entirely 
from the p T dependence of the Sivers function, and was further integrated out building the fit in terms of specific 
moments of the function itself. The x dependence of that approach is retained, but a different flavor-dependent 
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normalization and an explicit p T dep endence are introduced that are bound to the shape of the recent RHIC data 
on pp^ — > ttX at y/s = 200 GeV [2|j, where large persisting asymmetries are found that could be partly due to the 
leading-twist Sivers mechanism. The expression adopted is 

ftJ{x,p T ) = N f x(l-x) M2PlP J fl(x,p T ) 



2 



4 
2, 



where po — 2 GeV/c, and N u = —Nd = 0.7. The sign, positive for u quarks and negative for the d ones, already takes 
into account the predicted sign change of f^ T from Drell-Yan to SIDIS. 

Again, following the steps described in Sec. III-2 of Ref. [2l|, we can directly insert Eq. (|14() into Eq. and get 

( 2 Po q T \ 2 8N U + N d 
c<* Ml- 9 ■ ( 15 ) 

The q T shape is different from Eq. (|13fl and the peak position is shifted at larger values. This is in agreement with a 
similar analysis of the azimuthal asymmetry of the unpolarized Drell-Yan data (the violation of the Lam- Tung sum 
rule But, more specifically, it is induced by the observed x F — q T correlation in the above mentioned RHIC 

data for pp^ — > irX, when it is assumed that the SSA is entirely due to the Sivers mechanism. This suggests that the 
maximum asymmetry is reached in the upper valence region such that x F ~ x 2 ~ (<7t)/5 |2q. 



B. The Boer-Mulders effect 



Contrary to the Sivers effect, the lack of data for the Boer-Mulders effect does not allow to build reasonable 
parametrizations either of hf ^ (x,p T ) or of h{(x,p T ). Therefore, similarly to what was done in our previous papers [TtI 
|30j , the strategy of the numerical simulation is based on making guesses for the input x and p T dependence of the 
parton distributions, and on trying to determine the minimum number of events required to discriminate various SSA 
produced by very different input guesses. In fact, this would be equivalent to state that in this case some analytic 
information on the structure of these T MD p arton distributions could be extracted from the SSA measurement. 

Following the steps in Sec. IV C of Ref. [lj| and in Sec. VI of Ref. [lj|, the p T dependence of the parton distributions 
is parametrized as 



f{(x,p T ) = ^ e - a -rif{(x) 

7T 



h{(x,p T ) = ^ e - aT ^h{(x) , (16) 

7T 



where a T = 1 GeV -2 and M c = 2.3 GeV. In particular, the p T dependence of is fitted to the measured cos2</> 
asymmetry of the corresponding unpolarized Drell-Yan cross section, which is small for 1 < q T < 3 GeV/c (see, 
for example, Fig. 4 in Ref. 15]). Correspondingly, the sin(0 + ips) SSA will turn out to be small for the considered 
statistically relevant q T range (see Sec. IIIIBTi . 

Inserting the expressions (fTfi)l into Eq. (fTT|) . we get 



c 4 = 



2M C q T E/ 4 flixrJ/H^ h[{x 2 - f/H\) + (f «-» /) 

(17) 



q\ +4M2 J2 f e 2 f fl(x 1 ;f/H 1 )f((x 2] f/H 2 ) + (/ ~ /) 
2M c q T f(x 1 ;(f)/H 1 )h 1 (x 2 ;(f)/Hl) _ 2M C , q T hx(x 2 ; (f) / ' h\) 



ql + 4M 2 f( Xl ; U)/Hi) h(x 2 ; (f)/H 2 ) ~ <& + 4M 2 ] x {x 2 - (f)/H 2 ) ' 

where the second step is justified by assuming that the contribution of each flavor can be approximated by a corre- 
sponding average function [l7j . 

Two choices with opposite features will be selected for the ratio hi(x 2 ; {f) /H\ )/fi(x 2 ; (f)/H 2 ), namely the ascend- 
ing function ^/x~2 and the descending one \fl — x 2 , that both respect the Soffer bound. The goal is to determine 
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the minimum number of events (compatible with the kinematical setup and cuts) required to produce azimuthal 
asymmetries that can be clearly distinguished like the corresponding originating distributions. We identify this as 
the criterion to establish when information on the analytical structure of the involved parton distributions can be 
extracted from SSA data. 



III. RESULTS OF THE MONTE CARLO SIMULATIONS 

In this Section, we present results for Monte Carlo simulations of both the Sivers and the Boer-Mulders effects in 
the Drell-Yan process 7r ± pT _> ^+ ^~ X using input from the previous Sec. Ill Al and III 51 respectively. The goal is 
twofold. On one side, to explore the sensitivity of the simulated asymmetry to the different input parametrizations of 
Eqs. lfT2f and fl^jl . as well as to directly verify, within the reached statistical accuracy, the predicted sign change of 
the Sivers function between SIDIS and Drell-Yan [2(j. On the other side, to make realistic estimates of the minimum 
number of events required to extract as detailed information as possible on the chiral-odd distributions hj~ and hi. 

We consider pion beams with energy of 100 GeV hitting a transversely polarized proton target such that \fs ~ 14 
GeV, i.e. the same cm. energy available at HESR at GSI in the socalled asymmetric collider mode with antiprotons of 
15 GeV and protons of 3.3 GeV 01 ■ The transversely polarized proton target is obtained from a NH3 molecule where 
each H nucleus is fully transversely polarized and the number of "polarized" collisions is 25% of the total number of 
collisions |17j . The muon pair invariant mass is constrained in the range 4 < M < 9 GeV, in order to avoid overlaps 
with the resonance regions of the cc and bb quarkonium systems. At the same time, the theoretical analysis based on 
the leading-twist cross section (QJ should be well established, since higher-twist effects can be classified according to 
powers of M p /M, where M p is the proton mass. 

In the Monte Carlo, the events are sorted according to the cross section (0J, supplemented by Eqs. J5J) and ©. The 
asymmetry is simulated by Eq. (Q. In particular, for the Sivers effect we use Eqs. iJHJ and (I13f) or (|15H . according to 
the input parametrization selected for the Sivers function. For the Boer-Mulders effect, we use Eqs. I|10f) and i|17[l . 
The events are divided in two groups, one for positive values (U) of sin(0 — cj)g 2 ) in Eq. (JSJ or of sin(<^> + 4>s 2 ) m 
Eq. 1)11) |). and another one for negative values {D), and taking the ratio (U — D)/(U + D). Data are accumulated 
only in the X2 bins of the polarized proton, i.e. they are summed over in the x\ bins for the pion, in the transverse 
momentum q T of the muon pair and in their zenithal orientation 9. 

Proper cuts are applied to the q T distribution according to the different inputs. As for the Sivers effect, the flavor- 
independent Lorentzian shape in the p T dependence of Eq. I|12[l produces a maximum asymmetry for q T ~ 1 GeV / c 
and a rapid decrease for larger values. Consequently, transverse momenta are selected in the range 0.5 < q T < 2.5 
GeV/c, because for larger cutoffs the asymmetry is diluted. For the case of Eq. I|14fl . the peak position in q T is shifted 
at higher values and the cut is modified as 1 < q T < 3 GeV/c. In this way, the ratio between the absolute sizes of 
the asymmetry and the statistical errors is optimized for each choice, while the resulting (q T ) ~ 1.8 GeV/c is in fair 
agreement with the one experimentally explored at RHIC 26] . As for the Boer-Mulders effect, we keep the latter cut 
1 < q T < 3 GeV/c. The 9 angular dependence for the Boer-Mulders effect is constrained in the range 60° < 9 < 120° 
due to Eq. ijTU|) . because outside these limits the azimuthal asymmetry is too small 0. On the contrary, for the 
Sivers effect there is no need to introduce cuts because of the (1 + cos 2 9) term in Eq. JSJ plj . 

We have considered different initial samples. The Sivers mechanism is explored starting from 100 000 events with 
the ir~ beam and 25 000 with the tt + beam, because the Monte Carlo indicates that the cross section involving tt + 
is statistically disfavoured by approximately the factor 1/4 |37l ]: in such a way, the two samples can be collected in 
the same time. As for the Boer-Mulders effect, the lacking of any parametrization makes it impossible to perform an 
isospin analysis; hence, we used 50 000 events with the ir~ beam. Statistical errors for (U — D)/(U + D) are obtained 
by making 10 independent repetitions of the simulation for each individual case, and then calculating for each X2 bin 
the average asymmetry value and the variance. We checked that 10 repetitions are a reasonable threshold to have 
stable numbers, since the results do not change significantly when increasing the number of repetitions beyond 6. 

A. The Sivers effect 

In Fig. [21 the left panel a) displays the sample of 100 000 Drell-Yan events for the ir~p' i — > /i + /i~Y reaction at 
y/s ~ 14 GeV as they are collected in X2 bins for muon invariant mass in the 4 < M < 9 GeV range. The right panel 
b) contains 25 000 events for the tt + p^ — > fi + n~X reaction in the same kinematic conditions. Both samples can be 
accumulated approximately in the same time according to Eq. <|13[) based on the parametrization (|12|l of the Sivers 
function |25| : as already discussed, the transverse momentum distribution is constrained in the range 0.5 < q T < 2.5 
GeV/c. For each bin two groups of events are stored, one corresponding to positive values of sin(</> — </>s 2 ) in Eq. (jSJ) 
(represented by the darker histogram), and one for negative values (superimposed lighter histogram). Since the 
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FIG. 2: The samples of Drell-Yan events for the Sivers effect in the ir ± p^ — + /i + ff~ X reaction at y^i ~ 14 GeV, 4 < M < 9 
GeV, and 0.5 < q T < 2.5 GeV/c, using the parametrization of Eq. (II lit (see text), a) left panel: 100 000 events with the n~ 
beam; the darker histogram collects events with positive sin(0 — 4>s 2 ), the superimposed lighter histogram collects the negative 
ones, b) right panel: the same for 25 000 events with the 7r + beam. 

qq — ► 7* mechanism tends to populate the phase space for the lowest possible r values 0, 0, |3(J compatible with 
the explored range 0.08 < r = x\ x-2 < 0.4, this reflects in a x\— integrated distribution which is peaked for xi values 
in the valence domain. 

asymmetry 
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X 2 

FIG. 3: The asymmetry (U — D)/(U + D) corresponding to the histograms of Fig. El where U identifies the darker histograms 
and D the superimposed lighter ones (see text). Triangles for the parametrization of Eq. Ill 2ti using the n~ beam and with 
N u > 0; squares for N u < 0. Open triangles using the -k + beam and with N u > 0; open squares for N u < 0. 

In Fig. |3 the asymmetry (U — D)/(U + D) is shown for each bin x 2 between the events of the previous figure 
accumulated for the positive (U) and negative (D) values of sin(^> — (j>s 2 ) in Eq. (JHJ. Average asymmetries and 
(statistical) error bars are obtained by 10 independent repetitions of the simulation. Boundary values of x-i beyond 
0.7 are excluded because of very low statistics. The triangles indicate the results with the 7r~ beam obtained by 
Eq. I|l 311 assuming that f^ T changes sign from the parametrization 1)12(1 of the SIDIS data to the considered Drell- 
Yan [23| . For sake of comparison, the squares illustrate the opposite results that one would obtain by ignoring 
such prediction. Finally, the open triangles and open squares refer to the same situation, respectively, but for the 
7r + beam. The sensitivity of the parameters in Tab. [I] to the HERMES results for the Sivers effect, reflects in a 
more important relative weight of the d quark over the u one in the valence x^ range, with opposite signs for the 
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corresponding normalization Nf, f — u,d. Consequently, in the valence picture of the (jr~)ir + — p collision where 
the (uu) dd annihilation dominates, the SSA for the Drell-Yan process induced by ir + has opposite sign with respect 
to 7r~. Moreover, it has an absolute bigger size because the dd annihilations are weighted more than the uu ones. 
Apart for very low x-i values where the parton picture leading to Eq. becomes questionable, the error bars are 
very small and allow for a clean reconstruction of the asymmetry shape and, more importantly, for a conclusive test 
of the predicted sign change in fi T . 
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FIG. 4: The same situation with the same notations as in Fig. [5] but for the parametrization of Eq. 114H with 1 < q T < 3 
GeV/c (see text). 



In Fig. 01 the Drell-Yan events are shown in the same conditions and notations as in Fig. |3 i.e. in the left panel 

a) 100 000 events for the ir~p^ — > /i + n~ X reaction at y/s ~ 14 GeV and for 4 < M < 9 GeV, and in the right panel 

b) 25 000 events for the ir + p^ — > /i + fj,~ X reaction in the same kinematic conditions. The difference is that the events 
are now collected according to Eq. (|15fl based on the parametrization (|14fl of the Sivers function (21J; the cut in the 
transverse momentum distribution is now 1 < q T < 3 GeV/c. Again, the darker histogram refers to events with 
positive sin(c/) — <fis 2 ) m Eq. ©, while the superimposed lighter histogram to the negative ones. Similarly, the density 
of events is peaked for values in the valence domain because of the dominance of the low r portion of the phase 
space. 



asymmetry 
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FIG. 5: The same situation with the same notations as in Fig.|3 but for the parametrization of Eq. I|14fl (see text). 



In Fig. [21 the asymmetry (U — D)/(U + D) is shown for each bin X2 between the events of Fig. ^accumulated for 
the positive (U) and negative (D) values of sin(0 — 0s 2 ) in Eq. ((HJ. Notations are as in Fig. 01 the triangles indicate 
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the results with the 7r~ beam obtained by Eq. (|15fl using a positive normalization N u , which already accounts for the 
sign change of f^ T from SIDIS to Drell-Yan; the squares illustrate the results obtained by ignoring such prescription, 
while the open triangles and open squares refer to the same situation, respectively, but for the ir + beam. Again, 
the opposite normalizations of the two flavors u, d, determine the opposite SSA between the tt~ and the 7r + beams. 
But now in Eq. i|14f) the relative weight of u and d distributions is the same, hence the absolute sizes of the SSA 
are approximately the same irrespectively of the charge of the 7r beam. As already anticipated in Sec. Ill Al the q T 
distribution induced by the parametrization Ijl4|l is also related to the observed x F — q T correlation in the RHIC data 
for pp^ — > irX 26] , when it is assumed that the SSA is entirely due to the Sivers mechanism. This suggests that the 
maximum asymmetry is reached in the upper valence region such that x F ~ X2 ~ (q T ) /5 ~ 0.4 for the considered cut 
in q T , as it is confirmed in Fig. [S] Similarly to the case of the other parametrization, the statistical error bars are 
very small and allow for a detailed analysis of the (universality) properties of fi T . 
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FIG. 6: The sample of 50 000 Drell-Yan events for the Boer-Mulders effect in the 7r~p T — > /i + ^jT X reaction at y/s ~ 14 GeV, 
4 < M < 9 GeV, and 1< q T < 3 GeV/c (see text), a) left panel for the choice h 1 (x 2 , (/}/#] )//i(a;2, (f)/H 2 ) = y/1 - x 2 «/} 
represents a common average term that replaces each contribution in the flavor sum, for further details see text); the darker 
histogram collects events with positive sin((/> + 4>s 2 ), the superimposed lighter histogram collects the negative ones, b) right 
panel: the same for hi(x 2 , (/) / H^)/ fi(x 2 , (f)/H 2 ) = y/x2. 



B. The Boer-Mulders effect 

In Fig.|SJ a sample of 50 000 Drell-Yan events for the 7r~p T — > fx + fi~ X reaction at \fs ~ 14 GeV is displayed in X2 
bins for muon invariant mass in the 4 < M < 9 GeV range and for 1 < q T < 3 GeV/c. Events are produced by the 
Boer-Mulders effect contained in Eq. 1171) . where the left panel a) refers to the choice h±(x2, {f) / H\ )/fi(x2, (/) jH-i) — 
\J\ — X2 and the right panel b) to the hi(x2, (/) /H\)l fi(x2, {f) / H2) — \/x~2 one. Here, (/) means that each term 
contributing to the sum upon flavors is replaced by a common flavor-averaged parton distribution. Following previous 
notations, for each bin the darker histogram represents events with positive values of sin(</> + 4>s 2 ) m Eq. i|10|) an d f ne 
superimposed lighter histogram indicates the ones with negative values. Similarly, the density of events is peaked for 
X2 values in the valence domain because of the dominance of the low t portion of the phase space. 

In Fig. [7| the asymmetry (U — D)/(U + D) is shown for each bin X2 between the events of Fig. accumulated 
for the positive (U) and negative (D) values of sin(0 + <fis 2 ) m Eq. (|10f> . Triangles correspond to the y/l — X2 input 
function, open triangles to the y/x~2 one. Both choices respect the Soffer bound between h\ and /1 and have an overall 
normalization 2/3, which seems a reasonable expectation on the basis of lattice results and first SIDIS experimental 
data 0| . The error bars represent statistical errors only. As it is evident in the figure, the open triangles describe a 
SSA which statistically reflects the ascending trend of the input function yjxi , while it is not the case for the other 
choice. Despite the small error bars, which allow to state that both SSA are nonvanishing and to distinguish the 
two cases in the narrow range 0.1 < X2 < 0.3, it is not possible to conclude that specific information on the analytic 
dependence of h\(x) can be extracted from such simulation, contrary to what is claimed in Ref. [13. The reached 
statistical accuracy indicates that the size of the sample may not be responsible for such failure. Rather, the p T 




11 



asymmetry 

0.2 
0.15 

0.1 
0.05 





-0.05 

-0.1 
-0.15 

-0.2 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 

X 2 

FIG. 7: The asymmetry (U — D)/(U + D) corresponding to the histograms of Fig. El where U identifies the darker histograms 
and D the superimposed lighter ones (see text). Triangles for hi(x2, (f)/H\ )/fi (x2, (f) / H2) = \A — X2- Open triangles for 
h!(x2, (f)/Hl)/Mx2, (J) /Ha) = y/xi. 

dependence of in Eq. (| 1 fcil) induces the overall small size of the displayed SSA. Moreover, the Soffer bound and the 
flavor independence of the analysis further reduce the selectivity power of the final Monte Carlo output. In particular, 
the latter issue calls for specific parametrizations of h{(x,p T ) and hf * (x, p T ), whose unavailability reflects the lacking 
of experimental data for sin(^) + </>s 2 ) asymmetries in single-polarized Drell-Yan processes. 

IV. CONCLUSIONS 

In a series of previous papers 0, |3(j , we investigated the spin structure of the proton using numerical simulations 
of single- and double-polarized Drell-Yan Single-Spin Asymmetries (SSA) for the p^p^ — > p + p~X process as well 
as for the pp^ — > p + p~X one [2l|. We selected muon pair invariant masses in the range 4 < M < 9 GeV (and also 
12 < M < 40 GeV for the case of proton beams), where there is no overlap with the resonance regions of the cc 
and bb quarkonium systems and the elementary annihilation can be safely assumed to proceed through the qq — > 7* 
mechanism. In particular, the Monte Carlo was based on the Drell-Yan leading-twist cross section, because higher 
twists may be suppressed as M p /M, where M p is the proton mass. 

As for single-polarized reactions, two interesting contributions generate azimuthal asymmetries of the kind sin(</> + 
4>s) and sin((/> — (f>s), where <f> and 4>s are the azimuthal orientations of the plane containing the final muon pair and of 
the proton polarization, respectively, with respect to the reaction plane. The first one involves the convolution of the 
transversity hi, the missing piece necessary to complete the knowledge of the nucleon spin structure at leading twist, 
and the Boer-Mulders hi , another chiral-odd parton density which is most likely responsible for the violation of the 
Lam- Tung sum rule, the long-standing problem of an anomalous cos 2<f> asymmetry of the corresponding unpolarized 
Drell-Yan cross section [l^. The second convolution involves the socalled Sivers function f^p a "naive T-odd" 
partonic density that describes how the distribution of unpolarized quarks is distorted by the transverse polarization 
of the parent hadron. As such, fyp contains unsuppressed information on the orbital motion of hidden confined 
partons and on their spatial distribution inside the proton [l9| . 

In this paper, we have reconsidered the same scenario but for the n^p^ — ► p + p~X process at \fs ~ 14 GeV, that 
can be reached at COMPASS with pion beams of energy 100 GeV and transversely polarized proton fixed targets. As 
with antiproton beams, the elementary mechanism is dominated by the annihilation between valence partons (from 
p) and valence antipartons (from n). Taking advantage on the high statistics reachable with pions, in our Monte 
Carlo we have simulated both sin(</> ± SSA in the Drell-Yan cross section. For the Sivers effect we have used two 
parametrizations of f^ T : the one of Ref. [25|, which was deduced by fitting the recent HERMES data for the sin(0— cj)g) 
SSA in SIDIS 0h the one of Ref. [2lJ 1 which is constrained by the recent RHIC data for the pp^ — > 11X process at 
higher energy 26], when it is assumed that the SSA is driven by the Sivers mechanism only. The main difference 
is that the former displays an emphasized relative importance of the unfavoured d quark, and it gives an average 
transverse momentum (q T ) of the lepton pair lower than the latter. Consistently, we have built SSA by integrating 
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the q T distribution with adequate cuts, namely 0.5 < q T < 2.5 GeV/c for the former parametrization, and 1 < q T < 3 
GeV/c for the latter one. Results have been presented as binned in the parton momenta X2 of the polarized proton, 
i.e. by integrating also upon the antiparton partner momenta x\ and the zenithal muon pair distribution 9 with no 
further cuts. For the Boer-Mulders effect, since there is no such abundance of data and fits, we have used, as we 
did in Ref. |17|. very different input test functions and we have explored the sensitivity of the simulated sin(</> + cf>s) 
asymmetry within the reached statistical accuracy, integrating q T in the range 1 < q T < 3 GeV/c. Again, results have 
been presented as binned in X2 by integrating also upon x\ and 9, but with the further constraint 60° < 9 < 120° 
induced by the factor sin 2 9 which drives the angular distribution of muon pairs. 

Given the very different situations for the two analyses, also the goals are different. For the Sivers effect, the 
numerical simulation aims to establish the necessary statistical accuracy to distinguish different input parametrizations 
and to test the (universality) pro per ties of the Sivers function, in particular its predicted sign change when going from 
SIDIS to the Drcll-Yan process 20]. As for the Boer- Mulders effect, the goal is to make input guesses and to try 
to determine the minimum number of events required to discriminate various SSA produced by very different input 
guesses, that would allow to extract as detailed information as possible on the chiral-odd distributions h± and hi. 

In all cases, sorted events have been divided in two groups, corresponding to opposite azimuthal orientations of 
the muon pair with respect to the reaction plane (conventionally indicated with U and D), and the asymmetry 
(U — D)/(U + D) has been considered. Statistical errors have been obtained by making 10 independent repetitions of 
the simulation for each individual case and, then, calculating for each X2 bin the average asymmetry and the variance. 
For the Sivers effect, a starting sample of 100 000 events has been selected for the tt~ beam. Since, from the Monte 
Carlo, the cross section with tt + turns out statistically unfavoured by a factor 1/4 [3]}, we have reduced the sample 
to 25 000 events for the ir + beam in order to compare situations with the same " effective luminosity" . As for the 
Boer-Mulders effect, because of the unavailability of fits and isospin analyses, we have used 50 000 events with the 
7r~ beam. In all cases, the 1/t behaviour of the cross section, induced by the 7* propagator, has a twofold effect. It 
produces the highest density of events for bins in the valence domain, typically for x 2 ~ 0.3. At the given y/s, it also 
grants that the considered invariant mass range allows to explore the most populated portion of phase space, while 
avoiding overlaps with ranges where the elementary mechanism could be more complicated and the leading-twist 
analysis more questionable. The direct consequence is that, with a very large statistics of pions available, very small 
error bars are reached, except for boundary X2 values. 

The availability of different parametrizations of the Sivers function, obtained from independent sets of data, allows 
for a deep analysis of the flavor dependence of the resulting Drell-Yan SSA, as well as for a test of the universal 
properties of this parton density. It turns out that the asymmetry always changes sign when switching from the 7r _ 
to the 7r + beam, because in the valence picture of the (n~)ir + — p collision the (uu) dd annihilation dominates, and 
both the parametrizations here considered have weights with opposite signs for the u and d valence quarks. The 
parametrization of Ref. [25| , being deduced by SIDIS data for the Sivers effect , displays a more important relative 
weight of the d quark over the u, which reflects in a smaller absolute size of the SSA with the n~ beam with respect 
to the 7r + case. No such evidence is shown by the parametrization of Ref. [2]]], constrained by data for the pp^ — * nX 
process at y/s — 200 GeV (2^|, where also the higher (q T ) induces a maximum of the asymmetry at higher X2, typically 
X2 ~ 0.4. In both the considered cases, we have simulated the asymmetry assuming or neglecting the predicted sign 
change of the Sivers function when replacing the SIDIS with the Drell-Yan process 20]. The corresponding results 
have, of course, opposite signs, but, noticeably, the very small statistical error bars allow to clearly distinguish between 
one choice or the other extreme. We conclude that with the considered sample of events it is possible to perform such 
important test of nonperturbative QCD using pion beams and transversely polarized proton targets in the kinematic 
conditions that can be prepared at COMPASS. 

Unfortunately, for the Boer-Mulders effect the lack of data and parametrizations of the involved parton distributions 
forbids a thorough analysis. The p T dependence of h± is inherited by fitting the measured cos 20 asymmetry of 
the corresponding unpolarized Drell-Yan cross section; for the statistically relevant range 1 < q T < 3 GeV/c, the 
s'm((f> + <fis) asymmetry turns out to be small. We have further approximated the transversity distribution by using a 
"flavor-averaged" ratio between h\(x 2 ) itself and the unpolarized parton distribution fi(x2), and we have simulated it 
by integrating upon x\,q T , 9, and inserting very different input test functions of x 2 , one ascending and one descending, 
but all satisfying the general constraints (like the Soffer bound, that puts a strong upper bound on the size of hi). The 
small statistical errors allow to conclude that the resulting (U — D)/(U + D) asymmetries, though small, are certainly 
nonvanishing. But the displayed trends in X2 are very similar and do not reflect the very different inputs. Hence, we 
conclude that in the present stage a measurement of such an asymmetry would not help in extracting information on 
the transversity hi and/or the Boer-Mulders function . 
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